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ABSTRACT 

Initial results are presented to demonstrate the performance of a linear sequen- 
tial estimator (Kalman Filter) used to estimate a LANDSAT-1 spacecraft attitude 
time series defined for four scenes. Previous work has shown that a similar esti- 
mator valid for one scene could facilitate precision geometric correction (recti- 
fication) and multitemporal registration to pixel (or better) accuracies. With the 
revised estimator a GCP poor scene — a scene with no usable geodetic control 
points (GCPs) — can be rectified to higher accuracies than otherwise, based on 
the use of GCPs in adjacent scenes. Attitude estimation errors have been de- 
termined by the use of GCPs located in the GCP-poor test scene, but which are 
not used to update the Kalman filter. Initial results achieved to date indicate 
that errors of <500m (rms) can be attained for the GCP-poor scenes. Operational 
factors are related to various scenarios. 
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1. INTRODUCTION 


Since the launch of LANDSAT-1 , much effort has been devoted to the 
development of digital methods for the geometric and radiometric correction 
of multi spectral scanner (MSS) data collected by the LANDSAT series of space- 
craft, As early as 1973 (References 1-3), examples of MSS data corrected 
to sub pixel accuracies by all digital methods were generated and reported 
on in the literature. Subsequent advances included scene/scene registration 
to sub pixel accuracies of multi temporal data (References 4,5). 

A particularly successful technique employed to reduce the data pro- 
cessing load and yet achieve high accuracy entails the novel use of a linear 
sequential estimator (Kalman filter) for estimation of the satellite attitude 
state vector (Reference 6). Results have been previously reported regarding 
the evaluation by various methods of the estimator's accuracy for precision 
registration of full scene LANDSAT MSS data (Reference 4,5) and rectification 
of single scene data (References 1-3, 7), 

In the previously referenced work, a given scene was corrected using 
knowledge of the spacecraft position, attitude and the sensor scan mechanism. 

For a given site* the locations of features within the imagery, referred to as contro' 
points, were employed to refine estimates of the spacecraft attitude state vector 
defined for the same site. Registration was accomplished by means of essentially 
the same procedure, but with two small modifications: (a) measurement uncertainties 

were made much smaller than in the case of rectification (effectively giving less 
weight to a priori data); and .(b) control points were used even if their precise 
geodetic locations on the earth's surface were unknown. 

While it often is sufficient to employ a site specific estimator for 
rectification/registration purposes, the technique is dependent upon the 
availability of control points for the site in question. In some cases of 
interest, the subject site may have few or no usable geodetic control points 
(GCPs) the locations of which on the earth's surface are known to high accuracy. 
Examples of such cases include sites remote from civilization, sites containing 
featureless terrain, sites with poorly defined land/water interfaces, or sites 
which have not been surveyed to establish geodetic control. 

*A site refers to a framed MSS scene approximately 185 km on a side. 
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In this paper, we report on the results of a technique which can in- 
crease the geometric correction and absolute location accuracy of such GCP 
poor sites when nearby sites do have adequate geodetic control. The technique 
is based essentially upon an extended definition of the attitude time series 
of the satellite, to encompass N sequential scenes in the direction of satellite 
motion. 

It is not the purpose of this paper to present the details of the algorithm 
and design for the extended attitude estimator, per se. Rather it is the objective 
herein to present initial performance results of the estimator, and to identify 
certain features of the estimator's performance which may be of general interest. 
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2. PERFORMANCE EVALUATION METHOD 


Based on available data for the LANDSAT attitude control system and 
attitude measurement system (Reference 8), as well as experience derived 
from previous work with the single scene estimator (References 6 and 7), 
it was decided to configure an initial series of tests involving four 
sequential LANDSAT scenes. . Table 1 lists the scenes employed for this 
purpose, derived from path 16 of the Worldwide Reference System. There were 
20 geodetic control points (GCPs) available per each scene, and the GCPs 
were approximately uniformly distributed within each scene. 

The evaluation procedure proceded in two steps. First, for a given 
scenario, the m scenes (one, two or three in number) which were GCP poor 
were identified. The GCPs from the GCP poor scenes were not used to update 
the attitude estimator (they are, in fact, the statistical control points 
of Reference 7). Estimation errors were computed at the GCPs of the GCP 
poor sites, based upon updates of the estimator by means of GCP location 
data in the GCP rich sites (4-m in number). 

The second step of the evaluation procedure reported on herein entailed 
calculation of the estimation errors at the GCPs of the GCP poor sites based 
only on a priori SIAT data recorded on the Computer Compatible Tapes purchased 
from the EROS Data Center? The SIAT data includes attitude, altitude and sub- 
point position of the satellite for approximately 10 points in time spanning 
each 185Km x 185Km framed scene. 

A relative comparison of the estimation errors obtained from the two 
steps of evaluation was the basis for determination of the improvement in at- 
titude estimation accuracy achieved with the multiple scene attitude estimator, 
in contrast to performance achievable without the use of any GCPs from the site 
in question. This evaluation was undertaken for various scenarios of GCP rich/ 

GCP poor scenes. Table 2 lists the various scenarios evaluated herein. 

Upon commencing the evaluations, it was determined' that" the SIAT at- 
titude data for successive scenes was discontinuous. Figure 1 shows the roll data 
for a sequence of scenes obtained directly from a NASA-produced Bulk Image Anno- 

*S I AT data for a number of scenes is contained on a Bulk Image Annotation Tape 
(BIAT). The scene CCTs employed for the initial evaluations were obtained 
before SIAT data was recorded on CCTs. Thus, it was necessary in this study 
to employ the NASA BIAT tape data, which subsumes the currently provided SIAT 
data. 
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TABLE 1. LANDSAT-1 Scenes Employed for Multiscene Analysis 
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TABLE 2. Scenarios Used for Performance Evaluation 
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tation Tape (BIAT). It appears from a study of this and other such data that 
the SIAT attitude data is the result of a cubic time series fit to processed 
telemetry data for attitude, and that the fit is performed on a single scene 
basis; hence, there is a discontinuity in attitude in the overlap region of 
two framed adjacent along track MSS scenes. 

In order to initialize the multiple scene attitude estimator, a global 
cubic attitude time series was fit to the aforesaid attitude data obtained 
from the BIAT. Selected BIAT data in the overlap regions were disregarded in 
order to achieve a good least squares fit to what appeared to be the underlying 
attitude data*. Figure. 1 also shows the fit data in relation to the BIAT data. 

The resultant fit errors are shown in Table 3. 

GCPs were employed alternately from the GCP rich scenes, to update the 
estimator. That is, one GCP was employed from the first GCP rich scene, and 
then another from the second such scene, until the last GCP had been pro- 
cessed. This method was employed to avoid biasing the global time series on 
the basis of only one GCP rich scene in the sequence. With each GCP update the 
estimation errdrs at the statistical control points (GCPs in the GCP poor scenes 
which do not update the estimator) were evaluated. The errors at each statistical 
control point (SCP) were evaluated. The error at each statistical control point 
was computed as the along track, across track and altitude components of the vector 
difference in earth centered rotating coordinates of the known statistical con- 
trol point position and the estimated position for the same control point. The 
estimated control point position was computed on the basis of the previous at- 
titude state vector (resulting from the previous update of the state vector by 
GCP location data in one of the GCP rich scenes) and the location of the control 
point feature within the image. The control points are located to +1/8 pixel 
accuracy by a manual control point designation process. This process employs 4:1 
digital enlargement by cubic convolution resampling (Reference .9) of a control 
point feature displayed on an interactive image display system. 


* 


Processed attitude data, before polynomial model fitting, could not be obtained. 
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TABLE 3. Fit Errors over 4 Scenes 1, 2, 3, 4 
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3. RESULTS 


The results for two types of scenarios are presented. The first type 
of scenario corresponds to the case in which one or more GCP poor scenes 
lie between two or more GCP rich scenes. This interpolation scenario con- 
trasts to the extrapolation scenario, in which one or more GCP poor scenes 
precede or follow one or more GCP rich scenes. 

A. Interpolation Results 

Interpolation results are listed in Table 4 for the scenarios 
identified by scene in Table 1. The Bulk Error corresponds 
to rms errors (meters) measured at GCPs which do not update 
the estimator (SCPs), on the basis of the NASA BIAT attitude 
defined for the individual GCP poor scene of each scenario. 

A Priori Error is the rms estimation error (meters) measured 
at the SCPs in the GCP poor scene of each scenario, on the 
basis of the smooth global attitude time series (single cubic 
time series spanning 4 scenes). The remaining two columns 
are the a posteriori estimation errors for all the GCP rich 
scenes, and the GCP poor scene of each scenario. 

A total of 10 GCPs were employed, regardless of the number of 
GCP rich scenes. It was found that the estimator converged 
rapidly with no more than 10 GCP updates. 

Comparison of the Bulk Error data to the A Posteriori Error for 
the same GCP poor scene is a direct measure of the efficacy of the 
global attitude estimator in comparison with processing such scenes 
in the absence of GCPs (in said scenes). The scene errors are 
dramatically reduced in all cases. There is a small remaining 
altitude component of error, but along track errors are reduced to 
60% of the Bulk Error and across track errors are reduced typically 
to 50% of the corresponding Bulk Error values. 

B. Extrapolation Results 

Extrapolation results are presented in Table 5 for the scenarios 
identified by scene in Table 1. In this case, not all scenarios ex- 
hibit improved performance, and those with non-conti guous single GCP 
rich/poor. scenes are distinctly worse off than the corresponding Bulk 
Errors. Contiguous single GCP rich/poor scenes are improved somewhat 
and significant performance improvement is achieved with two or more 
contiguous GCP rich scenes contiguous with a GCP poor scene. 
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TABLE 4. Interpolation Results 
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Altitude 1 GCP Rich Scene 



TABLE 5. Extrapolation Results 
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AT = Along Track * Scene Not Utilized 

CT * Cross Track 0 GCP Poor Scene 

H ■ Altitude 1 GCP Rich Scene 



4. CONCLUSIONS 


In general, it was found that interpolation mode error statistics were 
stable, that is, not sensitive to GCP processing order or location error of the 
GCPs. On the other hand, the extrapolation mode tended to be relatively less 
stable if extrapolated from one GCP rich scene. Not unexpectedly, performance 
deteriorates with separation Of GCP rich and poor scenes in the extrapolation 
mode. However, extrapolation with two or more GCP rich scenes significantly 
increases stability and accuracy. 

Clearly, more work needs to be done to examine the effects of attitude 
modeling errors, number of scenes which can be processed and the distribution 
of GCPs, Nonetheless, a significant improvement in geometric accuracy performance 
for scenes without GCPs can be expected with the use of a global state vector up- 
dated by GCPs in nearby sites. 


123 




REFERENCES 


1. S. S. Rifman, "Digital Rectification of ERTS Multispectral Imagery," 

Symposium on Significant Results Obtained from the Earth Resources 
Technology Satellite-1, New Carrol ton, Md. , NASA SP-327, p. 1131 
(March 5-9, 1973). 

2. S. S. Rifman, "Evaluation of Digitally Corrected ERTS Imagery," Symposium 

on Management and Utilization of Remote Sensing Data, Sioux Falls, South Dakota, 
p. 284 (Obtober 29-November 1, 1973). 

3. J. E. Taber, "Evaluation of Digitally Corrected Images," Third Earth Resources 
Technology Satellite-1 Symposium, Washington, D.C., NASA SP-351 , p. 1837 
(December 10-14, 1973). 

4. R. H. Caron, "Evaluation of Full-Scene Registered ERTS MSS Imagery Using 
a Multi temporal /Multi spectral Bayes Supervised Classifier," Fourth Annual 
Remote Sensing of Earth Resources Conferences, University of Tennessee 
Space Institude, Tullahoma, Tennessee (March, 1975). 

5. S. S. Rifman, et al . , "Second Generation Digital Techniques for Processing 
LANDSAT MSS Data," Earth Resources Survey Symposium, Houston, Texas 

(June 8-12, 1975). 

6. R. H. Caron and K. W. Simon, "Attitude Time-Series Estimator for Rectification 
of Spaceborne Imagery," Journal of Spacecraft and Rockets 12, 27 (January 1975). 

7. S. S. Rifman and W. B. ATlendoerfer, "LANDSAT MSS Rectification and Registration 
Results Obtained by Using a Precision Attitude Estimator," paper present to 

The Flight Mechanics/Estimation Theory Symposium, NASA, Greenbelt, Maryland 
(October 29-30, 1975). 

8. LANDSAT Data Users Handbook, NASA Document No. 76SDS4258 (September 2, 1976). 

9. K. W. Simon, "Digital Image Reconstruction and Resampling for Geometric 
Manipulation," Symposium on Machine Processing of Remotely Sensed Data, West 
Layfayette, Indiana (June 3-5, 1975). 


125 



